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Abstract 
> , 

The relations between solutions of the three types of totally linear 
f--^ , partial differential equations of first order are presented. The approach 

^sO \ is based on factorization of a non-homogeneous first order differential 

^=^ ■ operator to products consisting of a scalar function, a homogeneous 

^^ . first order differential operator and the reciprocal of the scalar func- 

tion. The factorization procedure is utilized to show that all totally 
linear differential equations of first order can be transformed to each 
other, and in particular to a homogeneous one. 



1 Introduction 



The method for solution of Lagrange partial differential equation is well 
known, and is found almost in every text book on partial differential equations [| 
|3|, ffl. Our goal in this article is to show how the factorization of a non- 
homogeneous first order differential operator leads quite naturally to simple 
relations between the solutions of three related types of Lagrange equations. 
Let E be an open subset of 3? n , L be a continuous vector field in E, and 
denote by C k (E) the set of real valued functions which are continuously dif- 
ferentiable of order k on E. The continuous vector field L may be viewed as a 



differential operator |I] from C l (E) to C°(E). For each real valued continu- 
ous function q : E — > R there corresponds an operator L+q : C 1 (E) — > C°(E) 
defined by (L + q)ip = Lip + qip,\/ip G C l (E), where Lip is the Lie derivative 
of the function ip with respect to the field L, and qip is the usual product of 
two functions. Our goal in this work is to study the relations between the 
solutions of the partial differential equations 

(i) L<\> = [it) (L+q)iP = 0, (Hi) (L+q) X = b, 
where b is a continuous function on E. The approach followed here hinges 
on factorization of the first order non- homogeneous differential operator (L + 
q) to a product of a scalar function, a homogeneous differential operator, and 
the reciprocal scalar function. 

2 Factorization of a First Order 
Non- Homogeneous Operator 

Let 1] G C l (E) be a non-zero solution of the differential equation 

(L + q)r ] = 0. (1) 

Equivalently, r] is any element in the kernel of the linear operator L + q that 
is different from zero. As a first step we assume that rj has no zeros in E, and 
hence i]~ x exists and of class C 1 (E). The general case in which rj vanishes 
on a subset 5 C E will be considered in section 4. We start by proving a 
useful operator equality on which hinges the method of reducing one type of 
Lagrange equations to another. 

Theorem 1 In C X (E) the following operator equality holds 

rjLr)- 1 = L + q. (2) 

Proof: for every ip G ^(E). 

(r]Lr]~ 1 )ip = r](r]~ 1 L — r]^ 2 ( y Lr]))ip = (L — r]~ 1 (Lr]))ip = (L + q)ip. 
We have used equation (0) to make the last step. 



Corollary 1 Equality (W is equivalent to 

L = V -\L + q) V (3) 

which shows that all operators of the form (L + q) which are based on the 
same field L may be transformed to L, and accordingly to each other. 



Corollary 2 The equality (j^) shows that the left hand-side must not be de- 
pendent on the particular solution r\ of equation (j^), since its right hand-side 
is not. Therefore if £ is another solution of ^), then by equation t^) and 
a similar equation written for the solution £, we have ^~ 1 ii]Lri~ 1 ^ = L. This 
yields L^/rj) = 0. 



Corollary 3 From ^) we deduce that 

L k = r ] - 1 (L + q) k r ] and (L + q) k = nL^' 1 (4) 

where k is a non-negative integer. If L is invertible then the latter rela- 
tion holds for all integers. We assume in relation §Q) that L is a C k ~ 1 vector 
field, q is a C k ~ x function and n is aC k function. 

Corollary 4 // (^) holds then it is easily checked that (L + kq)r] k = 0, and 
hence 

L = r ] ' k (L + kq) V k (5) 

In general, and for any real number a, we have 

L =\ r) \- a (L + aq) | 77 \ a (6) 

Corollary 5 If Q is a real valued continuous function on E then 

r l -\L + q + Q)r l = L + Q. (7) 



Corollary 6 Take Q = —A (A G R) in corollary JI2) to obtain 

(L - A)^ A = & (L + q - A)(# A ) = (8) 

The last relation states that: if ip\ is an eigenfunction of the operator 
L belonging to the eigenvalue A then rjipx is an eigenfunction of the oper- 
ator L + q belonging to the same eigenvalue A. 



Example 1 Take L + q = £ + 2x : ^(R) -> C°(R). Since r] = e~ x2 is a 
solution of (Jfy, we have 

CI ~,2 CI _2 

— + 2x = e x —e x 

It is obvious that every complex number A is an eigenvalue of the operator 
4- to which an eigenfunction ip\ = e Xx belongs. In accordance with the last 
corollary, it is easily checked that A is an eigenvalue of the operator 4- + 2x 
to which the function iyip\ = e~ x +Xx belongs. 



3 On Totally Linear Partial Differential Equa- 
tions 

Let (xi, ...,x n ) be a global system of coordinates on the region E, in which 
L is expressed as 

n 

L=J2 a k(xi,---,x n )d/dx k , (9) 

fe=i 

where the components a^ are of class C° on E. We shall describe the partial 
differential equation 

(L + q) X = b (10) 

where q(x) and b(x) are continuous functions on E, as totally linear. The 
totally linear equation 

(L + g)^ = 0, (11) 

will be referred to as the non-homogeneous reduced equation corresponding 
to (0)> or simply as the non-homogeneous equation. Equation (|10D is a 
special type of Lagrange equation. The method of solution of Lagrange 



equation flT0|), and consequently equations ( |T0] ) and ( |TT] ) is well known [0]. 
However we aim here to utilize equality (||) to reduce the non-homogeneous 
equation (|TT|) to the homogeneous equation 

L4> = 0, (12) 

and express its general solution in terms of a particular solution and the 
general solution of (|12|) . Alternatively, to express the general solution of (|TT1) 
in terms of n particular solutions. The results we have just pointed to are 
expressed in the following facts in which we assume that rj is a solution of 
( |TTD on E and that it has no zeros on E. 

Fl. A function O is a solution on E of L<p = if and only if i/j = rj(j) is 
a solution of (L + q)ip = on E. 

The proof is a direct consequence of corollary 1 in the previous section. 

F2. Let Q : E — > 3? be continuous. By corollary 5 in the previous 
section, a function ipo is a solution of ( pT|) on i£ iff \l/o = ^V'o is a solution 
of (L + q + Q)\fr = on E. In a more familiar language to the subject of 
differential equation, the transformation \l/ = rjip reduces the last equation 
to®. 

F3. If rj and £ are solutions of (JTTj) then C,/rj is a solution of (T2). Indeed, 



from equation ( p]) which is satisfied by r^ and £ we get r]LC, = £,Lr], and hence 

L(£/r ] ) = r ] - 2 (r ] LZ-ZL V ) = 0. 

F4. The general solution of the reduced non-homogenous equation (|TT|) 
is given by 

"0 = V f(.<h.,—;<f>n-l), (13) 

where 

<Pi = <l>i(xi,—,Xn) (z = l,....,n-l) (14) 

are (n — 1) functionally independent solutions of (p~2|) . 

Proof. According to the standard method in solving Lagrange equation 
M , the general integral of the homogeneous equation ([12]) is given by 



= /(01,....,0 n _l), (15) 

where / is an arbitrary (^function in its arguments. Now if ip is a solution 
of (|ll|) then by F3 ip/rj is a solution of dl2|) , and hence it must be of the form 
(jI3[). It follows that the general solution of ( PD is given by (p~3|) . 



F5. If 

T](x!, ....,£„), rj j (x 1 ,....,x n ) (j = l,....,n- 1) (16) 

are functionally independent solutions of (|TT| ) then the ratios 

<f>j = Vj/V U = l,....,n- 1) (17) 



are solutions of (12) . It is easily seen that these ratios are functionally inde- 
pendent. Hence 

<f> = f(Vi/v,--,Vn-i/v) ( 18 ) 

is the general solution of the homogeneous equation (|T2|), and 



1> = Vf{Vl/V,-;Vn-l/v) ( 19 ) 

is the general solution of the non- homogeneous equation (|TTD . If i] n is a 
further solution of ( pTT| ) then by (|19|) 

Vn/V = fo(Vl/V, -;Vn-l/v), ( 20 ) 

where /o is some specified C l function. 

F6. Let xo be a solution of the totally linear equation (|I0|). If x is an Y 
other solution of ( |10|) then (L + <?)(x — Xo) — 0- Hence if) = x ~ Xo is a 
solution of (|TTD, and consequently must be of the form (0). It follows that 
the general solution of fllOD is given by 

X = XO + ^/(01,"-,0n-l)- (21) 

If a second particular solution xi °f (0) is given, then r] = xi ~ Xo is a 
solution of flTT|) , and therefore the general solution of C[T0f) can be written as 
follows 

(x ~ Xo)/(xi - Xo) = /(0i, •••., 0n-i). (22) 

If rj2 is a third solution of flT0|), then by (p2]), (% 2 — Xi)/(Xi ~~ Xo) is a solution 
of ( |12"D and has accordingly the form flT5l). If it is known (n + 1) solutions 
Xo,Xi, ,Xn of (|TDp, then 

7? = Xl-XO, 7/1 = X2 - Xo, ••••, r] n -i =Xn~X0 
are solutions of (0), and hence 

(23) 
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are solutions of ( O ) . If these ratios are functionally independent then 

X — XO _ ,-,X2^XXl_ Xn — XO N /rj.x 

Xi ~ Xo Xi ~ Xo ' Xi - Xo 

gives the general solution \ of (|10|) in terms of (n+ 1) particular solutions. It 
is clear of course that these (n + 1) solutions of ( ]TD| ) are certainly functionally 
dependent. If Xn+i is a further solution of ( |10D then by (0) (xn+i — Xo)/(Xi — 
Xo) is some function of the solutions (^3J) of equation (|TJ). 



4 The General Case 

We proceed here to study the general case in which the solution rj of equation 
([!]) is defined only on an open subset E r] C E, and vanishes on a closed 
subset 5 C E. This latter assumption embodies most interesting cases which 
are encountered in practical examples, such as rj vanishing at a finite or a 
countable set, or outside an open subset of its domain of definition. In the 
present case the operator r\Lr\~ x exists only on E v — 5. Hence the equality (Q) 
must be replaced by the inclusion relation (r]Lr]~ l C L + q); it is an equality 
only on E v — 5. 

Some minor modification has to be made so that the corollaries of theo- 
rem! in section 2 and the facts in section 3 conform to the new assumptions. 
As an example the fact Fl in section 3 must be modified as follows: 

Fl. Let rj be a solution of ([!]) on E v which vanishes on a closed subset 
5cE v . 

(i) is a solution of ( |T2"D on E$ =^ r/cj) is a solution of ( |TT| ) on E^ fl E v . 

(ii) -0 is a solution of ( |TTD on E^ =3- ip/r] is a solution of (|12j) on E^nE v — S. 

Proof: (i) on E^ fl E v , where rj is defined, we have 

(L + q){r]4>) = (j)(L + q)r] + rjLcj). 

If L(f) = on E^ then (L + q){r}(j)) = on E v n £^. 

(ii) On E^ H E v — 5, where ^/?y is defined and is of class C 1 , we have 
L(ijj/t]) = 0, as it was shown in F3. 

It must be noted that the last fact determines the smallest domains on 
which and ip are defined. These domains may be extended by continuity 
to larger domains. 



The remaining facts and corollaries can be modified in a similar way. 

Example. Consider the operator L = xd/dx + yd/dy + zd/dz which 
is continuous on 3ft 3 . The general solution of the equation Lip = is = 
f(y/x, z/x), where / is an arbitrary C l function. The equation (L+3/2)ip = 
admits the particular solution 77 = x -3 / 2 on E v = {(i,t/,z) eS 3 :i> 0}. 
The function = x/y is a solution of L0 = on E^ = {(x, y, z) e 3ft 3 : y 7^ 0}, 
and the function ip = i](j) = x~ l l 2 y~ l is a solution of (L + 3/2)ip = on 
{(x,y,z) G 3ft 3 : y 7^ 0,x > 0} = E v fl -E^. It is clear that we could have 
adopted rj = \x\~ , and (j)' — \x\ jy (since L is linear), to obtain a solution 
ip' = \x\~ ' y^ 1 defined on 3ft 3 — {x = 0, y = 0} = E v > fl E^ . If we replace in 
the given equation 3/2 by 1 , then 77 = 1/x, and hence 7/0 = 1/y is a solution 
of (L + l)i/j = on £ . 

Consider now the totally linear equation (L + l)x = 3x 2 . Four particular 
solutions of this equation are xo = % 2 , Xi — x<1 + a;_1 ) X2 = x 2 + y _1 , xs — 
x 2 + z~ l . It can be checked that the formula fl24]) yields, in accordance with 
(|2l|) , the general solution x = x 2 + x _1 f(y/x, z/x). 



5 An Algebraic Comment 

Let A denote the commutative real algebra formed by the set of all functions 
defined in E. The set of solutions of the homogeneous equation L<p = 0, 
namely KerL, is clearly a subalgebra of A. The set of solutions of non- 
homogeneous equation (|TTD is a coset of this subalgebra determined by a 
particular solution of (0). If 77 is a particular solution of fllTl) defined every 
where in E, then Ker(L + q) = 77 KerL. It is clear that KeriL + q) is a 
sub- vector space of A. The set of solutions of equation (pi]) is a coset of the 
sub- vector space KeriL + q) determined by any particular solution of flIU|). 
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